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Abstract. An overview of the most important experimental results obtained in 
the first three running years with the BRAHMS experiment from Brookhaven 
National Laboratory (USA) is presented. The design of the experiment permits 
to measure the interesting physical quantities in a large ranges of rapidity and 
transverse momentum. Therefore, properties of hadron production vs rapidity 
and transverse momenta are presente. 
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1. Introduction 

The study of the ultra-relativistic heavy ion collisions has as one of its goals the 
recreation, in laboratory, of the conditions after Big Bang, at a few fj,s. 

According to the present cosmological theories, at this moment in the early 
Universe, a deconfined state of quarks and gluons has been existed. To obtain such 
phase of the nuclear matter, in laboratory, a transition to the hadron to quarks 
and gluons is necessary. Therefore, high nuclear temperatures and densities are re- 
quested. The nucleus-nucleus collisions at relativistic and ultra-relativistic energies 
can offer such conditions. 

BRAHMS experiment was placed at two o'clock on the ring of the Relativistic 
Heavy Ion Collider from Brookhaven National Laboratory (USA) . This experiment 
was designed to gather information on the interesting physical quantities character- 
ising various emitted particles in ultrarelativistic nuclear collisions as functions of 
transverse momentum, px, and rapidity, y. 

The yields as function of rapidity offer information on the nuclear density into 
a given collision, as well as on the produced entropy. The information on the colli- 
sion dynamics and the thermalisation degree is obtained from the spectral shapes 
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of the interesting physical quantities and their dependencies on rapidity. Important 
information on the collision dynamics is obtained from collision centrality measure- 
ments, and the information from the times in the reaction is included in the high 
pT parts of the spectra. 

In the first three years of running different experimental data have been col- 
lected. In the present paper a few of the most important experimental results 
obtained in Au-Au collisions at y'sjvw = ISOGey and ,/snn = 200GeV, as well 
as in d-Au collisions at ^/SNN = 200GeV are included. The experimental data 
obtained in p-p collisions at ^/SNN = 200GeV are used as reference data. 

The experimental results include information on charged particle multiplic- 
ity, pseudorapidity and rapidity distributions, participants, energy density, nuclear 
temperatures and radial flow, antiparticle to particle ratios, Coulomb momentum, 
chemical potentials, entropy per baryon and high px suppression. All quantities 
offer information on the behaviour of the highly excited and dense nuclear mat- 
ter formed in the overlapping region of the two colliding nuclei, as well as on the 
possible formation of some new forms of the nuclear matter. Some comments and 
comparisons with different assumptions are included, too. 

Before the presentation of the experimental results, a short description of the 
BRAHMS experiment is given. 

2. The BRAHMS experiment 

The BRAHMS (B road RA ngc H adron M agnctic S pcctromctcrs) Experiment is 
one of the experiments at Relativistic Heavy Ion Collider (RHIC) from Brookhaven 
National Laboratory (BNL) [1-5]. 

In BRAHMS experiment two different regimes can be studied, namely: (i) a 
baryon poor region with a high energy density, created at mid-rapidity; (ii) a region 
at high rapidities, near the initial nuclei, very rich in baryons at relatively high 
temperature. The BRAHMS experiment well identified charged hadrons - tt"*", tt", 
K~^, K~ , p and p - over a wide range of rapidity and transverse momenta at several 
energies and beams available at RHIC can be measured. 

Taking into account the accelerator system structure and the size of the ex- 
perimental halls the BRAHMS spectrometers are small solid angle devices. They 
provide semi-inclusive measurements in very difii'erent experimental conditions. 

The experiment was designed with two moveable magnetic spectrometers, namely: 
(a) Forward Spectrometer (FS); (b) Mid Rapidity Spectrometer (MRS) [1-5]. The 
Forward Spectrometer covers the angular region 2.3*^ < 9 < 30^, and the Mid Ra- 
pidity Spectrometer covers the angular region 30° < 6* < 95°. The pseudorapidity 
ranges covered by the two spectrometers are: 1.3 < ry < 4.0, 0.1 < 77 < 1.3, respec- 
tively. To the two magnetic spectrometers three event characterisation detector 
systems have been added; they provide global information. The three detectors 
systems are the following: Beam-Beam Counters (BBC), multiplicity detector and 
Zero Degree Calorimeters (ZDC). Together, the three detector systems will provide 
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Fig. 1. A top view of the BRAHMS Experiment 



centrality coverage in the mid rapidity region, namely: —2.2 < i] < 2.2, in the 
region 3.2 <| 77 |< 4.3, as weh as at 0°. A top view of the BRRAHMS experiment 
is shown in Fig.l. 

3. Particle identification and momentum determination 

Particle identification can be done using different ways. First of all, the detectors 
from the two spectrometers can be used to combine the time-of-flight measurements 
in the two hodoscopes (HI and H2 in Fig.l) with measurements in the threshold 
Cherenkov counter (CI) and the Ring Imaging Cherenkov Counter (RICH). The 
40 scintillator slats in the first hodoscope and the 32 scintillator slats in the second 
hodoscope - instrumented with photomultiplier tubes (PMT) at each end - allow 
the particle identification in the momentum range 1-20 GeV/c, in Forward Spec- 
trometer. The time-of-flight array in Forward Spectrometer is place at 8.6 m from 
the nominal vertex position. Mid Rapidity Spectrometer uses for particle identifi- 
cation two time projection chambers (TPCl and TPC2), a dipole magnet (D5) and 
a time-of-flight wall (TOFW) with 83 scintillator slats, placed at 4.3 m from the 
nominal vertex position. In the two arms of the BRAHMS experiment the overall 
time resolution is around 90 ps. A second time-of-flight wall was introduced in 
2003. In the Forward Spectrometer the separation of kaons and protons is acheived 
in the momentum range p < 4.5 GeV/c, and in Mid Rapidity Spectrometer in the 
momentum range p < 2.2 GeV/c. The time-of-flight performances for MRS and for 
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Fig. 2. Timo-of-flight performance for MRS at 90 degrees (a) and for FS at 12 
degrees (b). Mass-squared spectrum in the FS (c). Mass spectrum in FS using 
the RICH 



FS are presented in Fig. 2. 

Another way is related to the momentum spectra and cfFectivc mass. For fixed 
rapidities the spectra for positive and negative particles have been obtained. 
The following relation has been used; 

m2=p2(^-l) (1) 

where p is the particle momentum, t is the time-of-flight and L is the flight 
distance. 

Particle momenta are determined by projecting the straight-line tracks as re- 
constructed in the two TPCs to the magnet and calculating the bending angle at 
matched tracks using an effective edge approximation. For obtaining momentum 
spectra the vertex dependent acceptance maps have been used. Reactions with a 
vertex in the range from -15 cm to +15 cm were used in the analysis. The momen- 
tum resolution obtained with the BRAHMS experiment ^ < 0.02 [2-5]. 

The acceptances of the two spectrometers for different identified particles is 
shown in Fig.3. It is important to stress here that the acceptances for the two 
spectrometers, at a given field, for positive charged particles, are equal to the ac- 
ceptances for negative charged particles at the opposite polarity of the given field. 
Therefore particle ratios can be obtained. 
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Fig. 3. Acceptance maps for BRAHMS experiment 



4. The pseudorapidity and rapidity distributions and charged 
particle multiplicity 

For obtaining global information the Beam-Beam Counters (BBC), multiplicity de- 
tector and Zero Degree Calorimeters (ZDC) have been used. The Beam-Beam 
Counters provide the initial trigger and vertex information. They provide the start 
time for the time-of-flight measurements, too. Each of the elements of the BBCs has 
3 cm and 4 cm UVT Cherenkov radiator read out by PMTs. The right side of the 
BBCs consists from 36 elements and the left side consists from 44 elements; they are 
positioned at ±2. 2m from the nominal vertex interaction. The multiplicity detector 
provides information on the collision centrality. 24 segmented Si-detectors with 168 
channels (SiMA = Si Multiplicity Array) and 40 scintillator squared tiles with 12 
cm length (TMA = Tiles Multiplicity Array) are included in this hybrid detector. 
The pseudorapidity range covered by the multiplicity detector is —2.0 < rj < 2.0. 
Common devices for all experiments at RHIC are the Zero Degree Calorimeters [6]. 
At the BRAHMS experiment they arc placed at ±14m from the nominal vertex 
interaction, situated behind the two DX beam-line magnets. For Au-Au collisions 
the two ZDCs offer information on the collision centrality measuring the forward 
going spectator neutrons. The coverage provides by the global detectors is in the 
region —2.2 < rj < 2.2, in the region 3.2 <| ?7 | 4.3, as well as at 0°. They have 
been used in the three runs to measure the overall charged particle production and 
to determine collision centrality. SiMA and TMA can be hit by many different 
particles, including background particles. Therefore, to obtain the multiplicity in- 
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Fig. 4. (a) Pseudorapidity distribution of charged particles in Au-Au collisions 
at 200 AGeV. (b) Multiplicity spectrum for Au-Au collisions at 130 AGeV 



formation it is necessary to convert the deposited energy into a number of hits. 
The energy loss signal is transformed to charged particle multiplicity dividing to- 
tal deposited energy to the expected deposited average energy by one particle in a 
tile. Conversion factors are evaluated from GEANT simulations. The background 
contributions are considered. They are between 10% for SiMA elements and 30% 
for TMA elements. All detector elements were calibrated in energy using pulsers, 
cosmic rays and sources measurements. 

The pseudorapidity distributions, j^, has been obtained from the measured 
energy depositions corrected using the above given corrections (Fig. 4a). Selecting 
appropriate ranges in the multiplicity spectrum (Fig. 4b), different cuts in collision 
centrality have been performed [7, 8] . Usually, the cuts in collision centrality are 
expressed in terms of the fraction of the nuclear reaction cross section by normal- 
ization to the integral of the TMA spectrum obtained, with a near minimum bias 
trigger. This trigger requires energy deposition in each of the two ZDCs above 25 
GeV. An additional condition is that at least one tile having a hit. Almost 97% of 
the nuclear interaction cross section can be selected by this requirement. 

In Au-Au collisions at ^^snn = 130 GeV the charged particle multiplicity on 
the whole pseudorapidity range was N^j^'^ = 4050 ± 300, and for the same collisions, 
at y^sWN = 200 GeV the charged particle multipHcity was A'^^^oo ^ 5100± 300. 

For the most central Au-Au collisions (0-5%), the values obtained from pseudo- 
rapidity distributions, at = 0.0, are: ^Ijf = 550±30, ^1^°° = 610±50, respec- 
tively. These values arc in good agreement with HIJING code predictions [10,9,4]. 
The full widths at the half maximum are: Aiyi^o = 7,5 ± o.7, Ar;^™ = 7.9 ± 1.0, 
respectively. 

It is important to stress that similar values are obtained in the full pseudorapid- 
ity range. This fact indicates the existence of a large plateau for particle production 
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Fig. 5. Invariant spectra of charged hadrons from Au-Au collisions at 200 AGeV 
for 2 pseudorapidities: (a) 0.0, (b) 2.2, respectively. The results for d-Au collisions 
are included, too. 



in very central collisions (6%). Also, this value for 77 = 0.0 is in good agreement with 
the value reported by PHOBOS Collaboration, value obtained in similar conditions 
[11]. 

The rapidity or pseudorapidity distribution can be used to estimate energy 
densities at the meson production, according Bjorken's formulae [12]: 

3 <pT> dN 

2 tttR^ dr) ^ ' 

For Au-Au collisions at ^JsWn = 130 and ^/snn = 200 the estimated energy 
densities are "i.&GeV/Fm^ and 4:.5GeV/Fm^, respectively. 

Taking into account that the pseudorapidity distribution resembles Bjorken 
boost invariant assumption hydrodynamic calculations can be considered further. 
Therefore, some considerations on momentum spectra and the inverse slope are 
included. 

5. Momentum spectra and some dynamic information 

Usually, the high momentum particles are considered as a good probe to investigate 
the conditions prevailing early in the evolution of a system. Such particles can be 
associated with different production mechanisms like jet production; they can loss 
energy due to induced gluon radiation passing through a medium with a high density 
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Fig. 6. Transverse mass spectra for charged pions, kaons, protons and antipro- 
tons for the most central Au-Au collisions at 200AGeV, at y=0 



of color charges. As a consequence, a depletion of the high transverse momentum 

component can by observed in the spectra. 

In this paper the positive identified hadron momentum spectra, for 2 pseudo- 
rapiditiies {rj = 0.0, r] = 2.2, respectively), as well as those of their antihadrons 
are included (Fig. 5). For each peseudorapidity different collision centralities are 
considered. These momentum spectra have been obtained as was presented pre- 
viously. The invariant inclusive spectra of negative hadrons and positive hadrons, 
respectively, at 77 = 0.0, permitted the estimation of the average values of the 
transverse momentum, for negative and positive hadrons. They are the following; 
< Pt >h-— OASOGeV/c, < pr >h+= 0A50GeV/c, respectively. These values are 
higher than those obtained by extrapolation at other energies [13,15]; for example, 
at ISR and UAl energies average transverse momentum is < pr >ex= 0.365Gey/c. 

5.1. Inverse slope versus mass and centrality 

The transverse momentum and transverse mass spectra have been used to estimate 
the inverse slope parameters or nuclear temperatures, T, other physical quantity in 
the study of the conditions in early in the evolution of a system. 

The transverse mass spectra for the hadrons detected with the BRAHMS ex- 
periment - for the most central Au-Au collisions at ^/snn = 200GeV(0-10%), at 
y=0 - are presented in Fig. 6. For fitting the invariant transverse mass spectrum the 
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Table 1. Experimental values of the temperatures for Au-Au collisions at 
200AGeV and y = 0.0 



following 



Centrality 


0-10% 


10-20% 


20-40% 


40-60% 


T(7r+) 


243.4±0.9 


245.5±1.1 


241.2±1.2 


230.3±2.0 


T(7r-) 


242.6±0.9 


241.6±1.1 


240.7±1.2 


231.5±1.9 


T{K+) 


303.6±3.3 


302.6±4.1 


292.5±4.5 


281.3±10.2 


T{K-) 


306.1±3.4 


310.6±4.2 


300.6±4.7 


288.5±9.0 


Tip) 


360.4±3.3 


353.6±3.6 


344.1±4.5 


317.4±9.8 




358.9±3.6 


342.0±3.8 


335.1±4.8 


306.2±11.4 


; relation has been used: 








1 


1 dN 


1 dN 


1 _ rriT — mo 


Nev 2'KmT drriTdy 


2ti dy T{T 


+ moY 


T 



(3) 



Here, ^ is considered as a free parameter, for different centrality cuts. The experi- 
mental results for the slope parameter T, for 4 different centrality cuts are included 
in Table I. The observed dependence of the tcnipcraturc on the particle mass (T 
increases with the increase of the hadron mass) , as well as the its dependence on 
the collision centrality (T increase with the increase of the collision centrality) is 
an indication of the radial expansion. Experimental values of the temperatures are 
higher than the values obtained from the simulations with HIJING and UrQMD 
codes [10,14]. 

From Table I it can observe that the pion's temperature is significantly higher 
than its rest mass. 

5.2. Collective transverse Row 

The mass dependence of the slope parameters provides evidence of the collective 
flow. The following relation between the slope parameter and particle mass has 
been proposed [15]: 

T = Tfo + m<v^> (4) 

where Tfo is the freeze-out temperature, considered as the temperature for that 
the particles cease to interact with each other, and < > is the average collec- 
tive flow velocity. The freeze-out temperature is related to the thermal, random 
processes, and the velocity is related to the collective processes. 

A fit with this formulae leads at freeze-out temperatures around 220 MeV, 
for all collision centralities, and collective flow velocities around 0.37c. Therefore, 
different relationships among these physical quantities have been proposed in the 
frame of different models [16, 17]. 
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Fig. 7. BRAHMS preliminary results for the most central collisions (0-10%) in 
comparison with the results from other experiments at lower energies 



Table 2. Freeze-out temperature and transverse collective flow velocity 



Ccntrality T/„ /3t 



- 


- 10% 


119±1 


0.626±0.005 


10 


-20% 


125±2 


0.602±0.006 


20 


-40% 


127±2 


0.583±0.008 


40 


-60% 


134±5 


0.525±0.023 



The experimental results for freeze-out temperature and transverse collective 
flow velocity, obtained in the framework of a blast wave model, are included in Table 
II. A very slow dependence of the freeze-out temperature on the collision centrality 
is observed (it decreases with the increase of the collision centrality) . The collective 
flow velocity has a more sensitive dependence on the collision centrality (it decreases 
with the decrease of the collision centrality). 

At RHIC energies the transverse collective flow velocity parameter is larger 
than that from collisions at lower energies (JINR Dubna, AGS-BNL, SPS-CERN). 
In the same time, the freeze-out temperature parameter is lower that obtained in 
nucleus-nucleus collisions at lower energies (see Fig. 7) [18]. 

6. Particle ratios 

6.1. General information 

Some information on the collision dynamics can be obtained from ratios of yields 
of antihadrons to hadrons, as well as from their dependencies on rapidity [19]. The 
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Fig. 9. Antihadron to hadron ratios as a function of rapidity. Errors bars 
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dependencies of the and | ratios on transverse momentum and collision 

centrality for Au-Au collisions at ^^snn = 200GeV - at 2 rapidities - are presented 
in Fig. 8. For a given rapidity no significant dependence of the ratios on transverse 
momentum and collision centrality is observed. Therefore, an analysis of the ratios 
dependence on the rapidity was done (Fig. 9). The different dynamical models do not 
describe such dependence. All ratios are essentially constant in the rapidity range 
y e [0, 1]. The values of /K^ and p/p ratios obtained in Au-Au collisions at the 
available energies at RHIC-BNL are highest values obtained. For example, at y = 
0.0, in Au-Au collisions at ^/SNN = 200GeV, for the most 10% central collisions, the 
following average values of the three ratios are obtained: < ^ >= 1.011 ± 0.006, 
< ^ >= 0.95±0.05and < | >= 0.75±0.04. Each ratio offers specific information. 

For example, the ^ ratio can be related to the Coulomb repulsion, the ratio can 
offers information on the strangeness content, and the | ratio can says something 
about the chemical equilibration (values of different chemical potentials). 

6.2. Coulomb interaction and the pion ratio 

The Coulomb interaction between net charge of the participant region and pions 
can be reflected by a modification in the transverse momentum. The Coulomb 
momentum is estimated using the following relation: 

, ^ idN"'' 1 

Pc^b-L - P-L,o =2e — 5 — — (5) 
ay Rf 

where p_L,o is the transverse momentum at freeze-out and p_L is the final momentum. 
The freeze-out radius is: 



Rf = Rgeom + 0.5-/3i • Tf (6) 

where Rgeom is the initial (geometrical) radius of the overlap zone at the collision 
time and (3 is the transverse flow velocity. 

Therefore, the pion ratio can be described using the following expression: 

where = yjiri"^ + {p±±pc)'^. 

The Coulomb effects in pion spectra are sensitive to the degree of stopping, 
the distribution of the positive charge and flow velocity of the participant region. 
The experimental results shown in Fig. 10 obtained in Au-Au collisions at reflect a 
reduced Coulomb effect {pc = 6.5 ± 1.9MeV/c). This behaviour can be related to 
the higher flow velocities of the nuclear matter for these collisions, in comparison 
with similar coUisions at lower energies {p^^^ = 20MeV/c, pf^^ = IQMeV/c). 
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Fig. 10. The pion ratio dependence on the transverse mass for the most central 
Au-Au colhsions at 200 AGeV (0-10%) - preliminary data 



6.3. Thermal interpretation of the kaon ratio dependence on the antiproton- 
proton ratio 

The baryon chemical potential can be related to the ratio ^ by the following relation; 

- = exp{ — ) (8) 

V T 

Using the experimental values of the ratio obtained in the BRAHMS experiment 
(0.64 ± O.OA{stat) ± 0.06(syst) for Au-Au collisions at ,/sWn = ISOGeV, 0.75 ± 
0.03(stat) ± 0.05{syst) for the same collisions at ^/snn — 20006^, respectively), 
as well as the following dependence of the baryon chemical potential on rapidity: 

+ (9) 

with a = 0.09 ± 0.01 and b = 1.69 ± 0.21, a increase of the baryon chemical 
potential with the increase of the rapidity has been observed (Table III) . The much 
reduced net-baryon density observed at mid-rapidity (10 ± 2) in comparison with 
forward rapidities (see Fig. 11) could explain this behaviour. Very different collision 
evolutions at different energies can be observed, also. 

For the ratio a similar dependence can be written: 

-— = e 3T T = (_) e T (10) 
K p 

In Fig. 12 the dependence: 



(11) 
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Fig. 11. Net-proton rapidity densities for central heavy ion collisions at AGS, 
SPS and RHIC 
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Fig. 13. The energy dependence of the chemical temperature 



is shown. For SPS experimental results the following results are obtained: 
a = 0.20 ± 0.01, and for BRAHMS experimental resuhs a = 0.24 ± 0.02. Both 
are obtained assuming T =170 MeV. This choice is in agreement with the results 
obtained supposing the following energy dependence of the chemical temperature: 

T,h - a - 6e~^^ (12) 

with a = 172.3 ± 2.8, b = 149.5 ± 5.7, c = 0.20 ± 0.03. For Au-Au collisions 
at y'SNN — 130Gey and ^/snn — 2006*6^ a common value is obtained, namely: 
Tch = 172.3±2.8MeT^. This results, together with experimental results obtained at 
lower energies (see Fig. 13), suggest that the chemical temperature saturates close 
to the critical temperature of 170 MeV extracted from lattice QCD calculations. 

It is important to stress here that the values of the baryon chemical poten- 
tial obtained in Au-Au collisions at the two energies (43.5 ± 3.2MeV, 28.8 ± 2.1 
respectively) are in agreement with the following dependence [21]: 

^s(s) = ^7j (13) 
1 + f 



Table 3. The dependence of the baryon chemical potential on the rapidity in 
Au-Au collisions at 200GeV 

y fiB[MeV] 

28.8±2.1 
1.5 46.7±8.9 

2.1 66.8±9.4 

3.1 123.1±10.9 
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with a = 0.967 ± 0.032Gey, b = 6.138 ± 0.399GeV. 

As general remarks on the dynamic information obtained from the behaviours of 
the antihadron to hadron ratios it is possible to affirm that the experimental values 
obtained in Au-Au collisions are the highest reported. The increase in ^ and ^ 
ratios can be consistent with the increase in collision transparency in comparison 
with lower energies. All ratios are constant, in the limit of the experimental errors, 
in the rapidity range y G [0,1]. This behaviour is in agreement with the large 
plateau observed in the charged particle rapidity distributions. The ratios permit 
to estimate baryon chemical potential and a chemical freeze-out temperature. At 
the RHIC top energy the value of the baryon chemical potential is lowest obtained 
up to now. At the same energy the baryon chemical potential increase significantly 
from mid-rapidity to forward rapidities. The net-baryon content is lowest at the 
mid-rapidity and could explain the behaviour of the baryon chemical potential. The 
correlation between and | ratios is well described by a statistical model. 

7. Some concluding remarks 

With the BRAHMS experiment from RHIC-BNL has been obtained a wealth of 
new and detailed information on relativistic heavy ion collision dynamics. These 
experimental data and results offer a tool to explore and quantify very hot and 
dense nuclear matter. The net-baryon density is very small (around 10), and the 
corresponding chemical potential is also small (around 29 MeV). 

The system exhibits a large transverse and longitudinal expansion with the 
azimuthally asymmetries being large, reflecting the initial partonic distributions. 
This behaviour is consistent with the assumption that the system has reached a 
hydrodynamic limit. This limit can be used to explore the equation of state of 
the hot and dense nuclear matter. Suppression of the high transverse momentum 
particles relative to elementary proton-proton collisions was observed in central 
Au-Au collisions. 

The heavy ion data from RHIC are consistent with formation of a hot dense sys- 
tem that exhibits hydrodynamic behaviour with rapid transverse and longitudinal 
expansion. Also, the system absorbs high transverse momentum probes correspond- 
ing to a large gluon density in the initially formed system. Taking into a account 
the value of the net-baryon density it is possible to affirm that is an almost baryon 
free system. 
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